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FOREWORD 


n.  .  Contract  F30602-71-C-0132,  the  Electro-Science  Laboratory  of 

t  ^ni^®rsity  has  been  Providing  direct  theoretical  support^ 
experimental  efforts  conducted  at  Rome  Air  Development  Center/  Both  the 
theoretical  efforts  of  03U  and  the  experimental  w^k  of  SdC  have  been 

p~iof:ri±8  the  ir^tion3  °r  57 

propagation  of  laser  signals  for  precision  pointing,  tracking  and 
paging  applications .  The  efforts  are  in  direct  rflpoSe  t^rwuire- 

Tl'ZoVvrZrT  AdVanCed  °Ptl0al  ^  “d  Tracking1^ 

al  theoretical  developnent  and  analyses  provided  by  03U  have  been 

aia,ed  at  providing  realistic  predictions  for  the  spatial  and  te^rora! 
behavior  of  propagated  laser  signals.  Analysis  of  averaging  time  re 

forr^Ud3oTOMrai?ted  th°  t0  detemtalll«  rlnimuLi  ttae  Requirements 
lor  vaiid  propagation  measurements  and  has  served  to  provide  systems 

d*“  ^  mll'fmizlne  operational  errors.  Ihe  extensioR  £  theRroSl 

P^auRR  ^RecRrf  U3SfUl  ^  thSt  U  provldes  »  measure  of  the 
cSlHo  obJaX  03363  r6  ^Petimehtal  results  are  dim- 
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ABSTRACT 


This  report  contains  detailed  information  on  the  theoretical  efforts 
of  Ohio  State  University  conducted  in  direct  support  of  the  RADC  Laser 
Propagation  Program.  The  efforts  described  are  concentrated  in  three 
main  areas:  the  examination  of  proper  averaging  times  required  for 
specific  propagation  experiments}  the  calculation  of  theoretical  curves 
predicting  results  of  various  experiments  along  slant  paths;  and  the 
calculation  of  predicted  results  for  the  temporal  spectra  of  particular 
atmospheric  effects. 

The  theoretical  development  and  calculations  have  been  aimed  at 
providing  useful  information  for  comparative  purposes  with  experimental 
efforts  and  for  use  by  systems  designers.  The  atmospheric  effects  in¬ 
cluded  in  the  analysis  are;  phase  structure  function,  angle  of  arrival 
and  coherence  length.  Both  the  spatial  and  temporal  characteristics 
have  been  considered. 
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INTRODUCTION 

This  is  Report  3163-3,  the  final  technical  report  under  Contract 
No.  F30602-71-C-0132  entitled  "Investigation  of  Laser  Propagation 
Phenomena."  The  effort  is  aimed  at  providing  theoretical  support  to 
the  RADC  Laser  Propagation  Program.  The  report  covers  the  period 
January  1,  1971  to  February  28,  1972. 

The  theoretical  support  is  in  the  area  of  linear  atmospheric 
propagation  phenomena.  Areas  of  interest  include  theoretical  support 
to  the  performance  of  propagation  experiments  and  the  interpretation 
and  processing  of  the  data  to  ensure  proper  match  between  theory  and 
experiment. 

The  work  presented  generally  relates  to  a  broader  program  aimed 
at  providing  basic  technical  information  on  atmospheric  imaging  and 
restoration  of  atmospherically  degraded  images,  on  pointing  laser 
beams  and  tracking  objects  viewed  through  a  turbulent  atmosphere. 
Experiments  are  in  progress  at  the  Rome  Air  Development  Center  to  provide 
the  desired  basic  information.  The  work  reported  here  first  provides 
theoretical  support  for  these  experiments  and  second  provide?  information 
useful  to  systems  engineers  who  might  use  the  information  generated  for 
particular  applications. 

This  report  describes  in  detail  work  performed  during  the  past  eight 
months  of  the  contract.  Work  performed  before  that  time  will  be  briefly 
summarized  subsequently. 

During  the  past  eight  months  work  has  been  concentrated  in  three 
main  areas:  the  examination  of  proper  averaging  times  required  for 
particular  propagation  experiments,  the  calculation  of  theoretical 
curves  predicting  results  of  various  experiments  along  slant  paths  and 
the  calculation  of  predicted  results  for  the  temporal  spectra  of  certain 
recorded  data.  These  areas  will  be  discussed.  There  will  also  be 
presented  a  writeup  on  a  suggested  redefinition  of  coherence  Tenths. 

The  work  on  averaging  times  is  in  direct  support  to  the  RADC 
experimental  program  because  it  supplied  procedures  for  the  per- 
formance  of  effective  experiments.  Indeed  data  has  been  examined  at 
RADC  amd  at  The  Ohio  State  University  to  check  the  recommended  procedures. 
Further,  an  averaging-time  meter  has  been  built  for  use  at  RADC. 

The  slant  path  curves  calculated  include  phase  structure  function 
and  angle  of  arrival  variance.  They  take  into  account  the  variation  with 
altitude  of  the  structure  constant  C^,  and  the  outer  scale,  L0.  These 
are  useful  for  systems  design  studies. 

The  temporal  spectrum  calculations  back  up  the  temporal  spectrum 
measurements  being  made  presently  at  RADC.  They  cover  the  time  lag 
spectra  of  the  phase  difference  temporal  correlation,  log  amplitude 
correlation,  and  phase  correlation;  and  'iclude  plane  and  spherical  waves. 
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Work  performed  during  the  first  portion  of  the  contract  is  covered 
in  reports  already  distributed  and  in  oral  papers  delivered  The  first 

Thatt^ennerannUalwteChni^al  report  (hereafter  designated  as  (3163-2) 

That  report  covered  a  review  of  averaging  time  studies,  a  biblioqraphv 

of  propagation  literature,  and  calculations  of  angle  of  arrival  correlation 
a5d  phase  structure  function  for  select  cases  whe?e  the  structure  ^rlmete? 
Cn,  and  the  outer  turbulence  scale  varied  with  altitude.  The  second  report’ 
was  a  study  of  angle  of  arrival  fluctuations.  It  presented  extensive 

correlations  a^fnrT^^  vfiation1for  a"9le  Sf  arrival  two  ooint 
correlations  and  for  large  aperture  angle  of  arrival  variance  In  it 

there  was  also  demonstrated  a  coincidental  identity  between  the  two  which 

greatly  simp  If led  calculations.  There  were  also  five  oral  papers  pre- 

sented,  one  (Hunt,  1971a)  at  the  Spring  Meeting  of  the  Optical  Society  of 

America  in  Tuscan,  Arizona,  April  5-8,  1972,  and  the  other  four  at  the 

0ctoberm5et1971Ol|whe  °f  America  in  Ottawa,  Canada, 

(ZintsmasieI-!l97l!yi971A):  ("yn9aard>  ,97,b>-  '97'b). 

contained  prTfc!int  ’TP°rt  is  dl'vided  into  three  separate  sections  each  self- 
contained.  The  sections  cover  averaqinq  times  slant  nath  raimiat-innr 
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I.  AVERAGING  TIMES 


data  InAshind?r^pSn  '^.Prese"t.work  on  averaging  times  for  propagation 
aata.  As  indicated,  this  work  is  aimed  at  supporting  the  RADC  exoeri 

mental  program  by  providing  means  for  assuringthat  a  given  data  Record 
is  sufficiently  long  so  that  assumptions  made  in  the  theoretical  de¬ 
velopments  are  born  out  experimentally.  Typical  criteria  are  that 
S  lit  ensemble  averages  be  identical  (the  ergodic  hypothesis) and 
what  the  average  is  independent  of  starting  time  (stationaritv)  The 

fourth1!!9  tlt6S  als°  depend  innately  on  the  quantity  being  measured 

d?fWnfm?!mLrt2UirHn9  °nger  times  than  second  moments,  and  requiring 
different  times  than  do  spectra,  etc.  This  work  is  essentially  an 

(3163-2)"  °f  that  reported  in  the  first  semiannual  technical  report 

,n  P?st  eight  months  averaging  time  work  has  progressed 

data  recordIdeartnRADfelhP  th®  f  am!natl'™  of  °Ptical  Phase  difference 
qive  a  real  Simp  construction  of  an  averaging  time  meter  to 

avpraninn  indi.cation  of  averaging  time  and  the  calculation  of 

be  SnsideS?5  ^  for  tefnporal  sPectra-  These  topics  will  now 


A.  Examination  of  RADC  Data 


Hiffprpn?o9U^V6’  19?I  SGVera1  runs  were  made  measuring  the  phase 
inference  at  two  points  on  an  atmospherically  degraded  spherical  wave 
Data  was  recorded  with  several  values  of  separation  of  the  two  do in?s 
The  data  was  computer  processed  at  RADC  to  remove  ambiguities  of  ' 

CopiesSofnthefln1rZed  *5  t0  9ive  means’  van'ances  and  temporal  spectra 
av£ragi^  times  d  ^  transmitted  to  osu  for  examination  of 

reDorthf3?63a?!S  anal^sl's  are  based  on  a  Previous  technical 

™  ,  6pi:  Ih  Procedure  was  applied  to  the  data  directly  and  to 

varianrp1^6  *he  data+*°  the  averaging  times  for  the  mean  and  the 

anance.  A  given  section  of  data  reduced  to  zero  mean  was  divided 

f^°eachaLbse«?onUbrCt™f  (desi9"ated  *>*  subscript,  s).  The  variance 
ror  eacn  suosection,  B  ,  was  computed  as  was  the  variance  for  the 

V  The  avera^e  of  the  subsection1  variances  and  the 

Sere  f,na9ny1Jalcru?atede.N’  Wee"  a"d  the  C™p,ete  sect1°"  fiance 


0) 


~N 


n  N  > 


N 


The  percentage  difference  was  then  plotted  as  a  function  nf  thp  ti™ 
duration  of  the  subsections.  A  typical  curSe  is  plotted  in 

tendpd^t  S!i-~CtT  duration  was  too  small  the  subsection  variances 
tended  to  differ  from  the  complete  section  variance.  However  as  the 
section  length  increased  the  subsection  variances  approached  more  dosely 
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the  S  ^at^ct^n*  "?  °",y  °"e  »■>«« o n, 

the  last  point  was  necessarily tlrl ^5"  variances  were  identical  so 
sections  of  increasing  duration  until  Ihn?nfr0CehUre  WuS  rePeated  with 
was  obtained  where  the  oercenJaop  0t.Sl?Ch  as  the  one  in  Fig.  1 

variance  and  the  complete  section  Jaria^p"^6^66"  the  mean  subsection 
sect, 'on  duration  mucE  «>r  sub- 


Fig.  1, 


vs.  averaging^t?men(30^sec^erecord)^nd  ensemble  avera9es 
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One  would  then  preselect  a  desired  degree  of  precision  and  look 
for  the  subsection  duration  giving  the  precision.  That  duration  would 
then  be  used  for  all  succeeding  runs. 

The  plot  in  Fig.  1  has  two  interpretations  (31 63-1 ) .  The  first  is 
that  the  averaging  time  is  sufficiently  long  that  the  variance  of  the 
data  has  come  within  a  given  fraction  of  the  stationary  value.  The 
second  is  that  the  fractional  difference  between  the  time  average  value 
and  the  ensemble  average  value  has  also  come  within  the  given  value. 


Figure  2  shows  a  plot  of  a  section  of  data  used,  redrawn  from  a 
computer  printout  of  the  digitized  values.  Each  point  represents  a 
tenth  second  average.  Figure  1  shows  a  typical  plot  of  the  percentage 
difference  ew  for  a  data  section  of  30  seconds  duration.  There  it  is 
seen  that  an  averaging  time  of  approximately  10  seconds  would  suffice 
for  a  value  of  eN  %  .03.  This  value  was  checked  by  increasing  the  data 
section  duration  up  to  as  much  as  five  minutes  with  the  same  result. 
Further  the  value  was  also  checked  by  delaying  the  starting  instant  of 
the  section  with  the  same  result.  The  conclusion  is  thac  for  the  mean 
value  alone  a  ten  second  section  of  data  would  have  been  sufficient  and 
that  the  variance  value  is  stationary  after  that  time. 

The  same  procedure  was  applied  to  a  second  data  section  formed  by 
squaring  the  original  phase  difference  data.  Figures  3  and  4  show 
typical  plots  of  the  percentage  difference  for  the  squared  data,  one 
for  twenty  and  the  other  for  one  hundred  and  twenty  seconds  section 
duration.  Both  show  the  expected  monotonically  decreasing  behavior. 
However  in  the  longer  section  the  normalized  difference  has  dropped  to 
0.0012  in  half  the  section  duration,  (60  seconds  point  on  time  scale), 
while  the  shorter  record  has  dropped  only  to  .04  at  half  section, 

(10  sec  point  of  time  scale)  duration.  Thus  it  appears  that  the 
averaging  time  of  sixty  seconds  would  be  sufficient  for  a  relative  un¬ 
certainty  of  .035  2:  v/.001 2 ,  a  value  well  within  the  five  minutes  of  data 
actually  recorded. 


B.  Averaging  Time  Meter 

It  is  certainly  desirable  to  have  all  data  records  of  sufficient 
duration.  On  the  other  hand,  excessively  long  data  records  give  in¬ 
efficiency.  Hence  it  would  be  desirable- to  have  some  method  of  de¬ 
termining  how  long  a  data  section  should  be  while  data  is  being  taken. 
Towards  this  end  two  schemes  have  been  proposed.  The  first  is  the  use 
of  an  on-line  digital  computer  to  perform  the  calculations  indicated  in 
the  previous  section.  The  second  scheme  is  an  electronic  technique  to 
provide  the  same  information  with  considerably  less  hardware  investment. 
In  its  essence  the  averaging-time  meter  computes  in  analog  fashion  the 
quantities  needed  to  determine  averaging  time  and  combines  them  so  as. 
to  read  out  on  a  meter  the  requisite  data  duration.  This  averaging  time 
meter  will  now  be  discussed. 
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Phas 


We  start  with  the  simple  equation  for  averaging  time  given  in  Eq.  (2) 
(3163-1),  Lumley,  (1964) 

(2,  T  =  §  ^1 

2 

Where  f  is  a  random  variable  and  e  is  the  mean  square  difference  between 
the  ensemble  and  time  averages  normalized  to  the  square  of  the  mean  value. 
I  is  the  integral  scale  of  the  autocorrelation. 


(3)  1  =  bToJ 


J  B(t)  dT 
0 


(4)  B(t)  =  (f(t)  -  f)(f(t+T)  -  7) 

The  restrictions  in  Eqs.  (3)  and  (4)  are  that  the  averaging  time  be  much 
longer  than  I,  and  that  the  data  be  normally  distributed.  Assuming 
a  gaussian  random  process,  then  the  averaging  time  is  (Lumley,  1964) 

>oo 

B(r)dT 

(5)  T  =  T-  0  w 

e 

In  terms  of  frequency  integrals  rather  than  time  domain  integrals 
Eq.  (5)  becomes 

(6)  T  =  - 

'  '  »oo 

e2  W(f)df 
J0 


Equation  (6)  contains  the  ratio  of  two  power  levels,  the  power  at  DC 
and  the  broadband  power.  T':us  the  problem  becomes  one  of  measuring 
these  two  quantities.  The  blcck  diagram  of  a  circuit  for  accomplishing 
this  is  shown  in  Fig.  5.  There  we  see  the  signal  entering  at  the  left 
and  split  into  two  parallel  paths.  The  top  path  squares  the  signal  and 
then  integrates  it  to  obtain  the  average  power.  This  implements  the 
denominator  of  Eq.  (6).  In  the  bottom  string  the  signal  is  integrated 
to  give  the  average  signal  which  is  then  squared  to  give  the  average 
DC  power  corresponding  to  the  numerator.  The  logarithms  of  the  two 
signals  are  then  subtracted  to  obtain  a  voltage  proportional  to  the  log 
of  the  averaging  time.  A  meter  with  a  log  scale  is  used  t.o  read  out 
the  averaging  time  directly. 

The  circuit  indicated  in  Fig.  5  has  been  constructed  and  is  shown 
in  Fig.  6  but  it  has  yet  to  be  tested.  Tests  are  expected  with  the  same 
random  signals  used  for  the  analysis  in  the  preceding  section. 
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C.  Temporal  Spectra  Averaging  Times 


It  is  of  interest  to  examine  the  temporal  spectra  of  random  data 
as  well  as  the  time  averages.  Indeed  Section  III  of  this  report  contains 
analytical  predictions  for  the  form  of  the  temporal  spectra  of  log-ampl itude, 
phase  and  nnase  difference  fluctuations  for  turbulence-degraded  light 
beams.  It  is  the  object  of  this  subsection  to  consider  one  averaging 
time  prediction  procedure  (Blackman,  1958)  as  applied  to  phase  difference 
temporal  spectra.  The  prediction  will  then  be  interpreted  to  show  that  it 
provides  one  limitation  on  averaging  times,  but  not  necessarily  the 
desired  limitation. 

The  discussion  is  based  on  the  expression  for  power  spectral  averaging 
times  found  in  the  literature  (3163-1),  (Blackman,  1958).  The  pertinent 
expression  is 


t;  -A- 

where 


(8)  We 


[ 


W(f)df]2 

fw2(f)df 


Tp  is  the  required  averaging  time,  e  is  the  relative  mean  square  un¬ 
certainty  in  the  power  spectral  measurement  and  We  is  a  mean  effective 
bandwidth.  Analytical  expressions  for  the  power  spectra  for  various 
propagation  variables  are  calculated  in  Section  III.  These  can  be  used 
to  obtain  analytical  estimates  for  the  mean  bandwidth  and  therefore  the 
averaging  time  using  Eqs.  (7)  and  (8).  We  present  such  a  calculation 
for  a  typical  case  of  interest,  phase  difference  measurement. 


For  phase  difference  measurements  the  power  spectrum  W*  (f)  is 
predicted  in  Section  III  to  have  the  form  6s 


(9a)  Wfis(f)  =  K  •  X-8/3 


where 

(9b)  K  =  .033  x  (2tt)8/  0  x  i 


" 

- 

n  sin  X 

"  X 

>♦(! 

■J 

•4/3 


(9c)  X  =  2irfr/v 


(9d)  R  =  1.077P/Lo 
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velocity  ^arun>  v  is  the  “*"«  S°mponent 

turbo, encc  structure  parameter,  and  ,  .  2,%  ?s wle^th. 

time  Ihe9expresslonP,inaEPr0?nafe  eStin,?te  °f  the  Pretlicted  averaging 
in  three  dSferen?  region^/  ’  "*  ,'ep,aCed  by  itS  ,Wti"9  *«!»«» 


(10)  K  =  X2/6R8/3  o  .<  X  <  R 

R  <•  X  ^ 

6  «  X 


=  X-2/3/6 
-  x-8/3 


express  iorTfo^w^is  given^  10  Eq‘  ^  W6re  performed-  The  resulti 


ng 


(11)  W  1  =  (2ttP/v) 


r  4r-173 

6-1/6 

— 

T3 

(gi 

4 

R,/.l)  2 

lire — 

9 

we  see  that^or^/L^l^th^ranqe^Ol 1<’  /L  ^"n  E,q,'/^  and  Fig-  7> 
range  1  <  WsV(P/v)  <  5  if IJ 1  "  p/L°  1  l’.0*  weV(p/v)  is  in  the 
and  wind  vefocity,  v  =  i5  mph  f%7  •>  m/O^nT^^n6^311^5  for  seParati 
the  range  of  values  for  the  decor^e^Hnn  I?  P  VV  m  and  e  =  0.1.  the, 

.33  seconds,  and  the  averaging  time  is  between’l3  Vand  fifitWGen  :°66  tI!' 
frequency  resolution  consistent  with  i,!  ?6'  i3,3  d  66  seconds-  Th< 
to  15.0  Hz.  consistent  with  the  desired  precision  goes  from  3.1 

may  noTtette'best’fo?  ^£.3«^'VaS^-0,\?ne  criterion  "Mel. 

the  decorrelation  time  and  noJ  Sn  deseed  soe^^,,y  \*  Is-  based 

cases  the  spectral  resolution  requirement  wif^hf  tK6S°  Utl°n*  In  many 
requiring  longer  data  duration  than  predict by Eq [  (7)^  $triCt  °ne 

section  of  data  T^econd^in^uration^JH  imfgi!?e  that  we  start  with  a 
spectrum.  It  will  have  a  resnlutinn^^^/oT  uU  ate.the  assoclated 
random  there  will  probably  be  scatter  in  thl  f2’  Since  the  data  is 
We  could  then  group  adjacent  points  Wtw  fr?quency  spectrum  points, 
to  reduce  the  scatter.'  ^his  inc?eas£9tL If£< J  aver5ge  within  a  group 
by  the  uncertainty  principle  that  is  mShiff  -1Vf  frequency  bandwidth 
data  section  into subsections f1nd1n^£1Jn?!1Vale!t  to  dividi"q  the 
and  averaging  the  subsecJiSn  spec?^?9  SP6CtrUm  °f  6ach  subsecti°" 
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Normalized  reciprocal  bandwidth  as  a  function  of  the  normalized  separation. 


The  interpretation  of  We  is  that  is  provides  a  limiting  largest 
bandwidth.  It  will  do  no  good  grouping  spectral  points  with  a  bandwidth 
larger  than  We  because  the  individual  points  would  be  no  longer  independent. 
That  is,  grouping  the  data  points  into  a  larger  frequency  interval  wojld 
be  equivalent  to  dividing  the  time  data  section  into  subsections  each 
shorter  than  the  correlation  time,  and  therefore  into  subsections  which 
are  no  longer  independent. 

If  better  frequency  resolution  is  desired  with  scatter  reduced,  then 
a  longer  section  of  data  must  be  used  to  provide  a  greater  density  of  points 
in  the  spectrum  and  hence  groups  of  smaller  bandwidth  but  the  same  scatter. 
As  indicated  this  is  more  apt  to  be  the  limiting  factor  on  averaging  time 
than  is  the  coherence  time  limit  indicated  in  Eq.  (7). 

As  an  example,  in  the  case  considered,  if  it  is  desired  examine 
frequencies  below  1.0  Hz,  than  a  0.1  Hz  resolution  would  be  more 
appropriate  than  the  range  3.0  to  15  Hz  cited.  One  must  go  to  a  data 

section  of  33  min.  duration  and  group  to  provide  more  spectral  Doints  to 
average.  K 

Summary 

To  summarize  the  results  of  this  section:  three  areas  associated 
with  averaging  time  have  been  considered.  In  the  first  area,  an 
examination  of  a  section  of  RADC  propagation  data  on  phase  difference 
fluctuations  indicated  that  sufficiently  long  averaging  times  had  been 
chosen.  In  the  second  area,  a  real  time  averaging  time  meter  had  been 
designed  and  constructed.  In  the  third  area  it  was  indicated  that  the 
limitation  on  averaging  time  for  spectral  data  was  probably  the  soectral 
resolution  desired  rather  than  the  percentage  uncertainty  required. 
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11  *  SLANT  PATH  calculations 

A.  Introduction 


amv; —< » --  .-sfjat-  ssvixisiZ' 

turbulent  1atmosphere^wMch^?ead<;tte  J^ufation^transfei^function^of^h 

“asaSas-iSa- 

curves  obtained  by  numerical  in^h  ,.-nd  Part  C  Presents  typical 

c°.fus?oen",CU,atii  -  -J  E™rS&  oTtte  l^g** 

The  Theory 


medium 


™^fsng?venaby  PhaSe  StrUCture  f°>-  *  locally  homogeneous  isotropic 
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where  p  is  the  separation  of  the  points  under  consideration 
k  is  source  wave  number  in  the  unperturbed  medium  =  ^ 

L  is  the  range 
<  is  the  spatial  frequency 

<J>n  is  the  index  of  refraction  structure  function. 

The  phase  structure  function  equation  is  simply  changed  to  the 
wave  structure  function  D  (p)  by  removing  the  cosine  squared  term  in 

Eq.  (1).  W 

In  the  case  of  the  larger  separations  the  cosine  term  does  not 
differ  from  unity  during  the  significant  part  of  the  integration, 
then  the  phase  structure  function  equals  the  wave  structure  function. 
For  small  separations,  the  cosine  term  oscillates  rapidly  during  the 
significant  part  of  the  integration  and  DS(P/L0)  is  1/2  DW(P/L0).  This 
is  identical  to  the  conclusion  for  the  plane  wave  case. 

The  elevation  angle  of  arrival  correlation  function  is  related 
to  the  phase  structure  function  by4 


(2a)  B  (p)  =  -!* 
“  2kd 


92Ds(p) 


and  the  azimuth  angle  correlation  function  by5 


(2b)  Bb(p)  = 


3DS(P) 


The  angle  of  arrival  variance  is 


(2c)  Bo(0)  =  Bp(0). 

The  Von  Karman  Index  spectrum6,7  including  inner  scale  i; 

/ ...  i.  \2 


<t>n(<»n,L)  =  .033  C^(n) 


<L  _"o 
n  5.92 


"  BMW 
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Cn(n)  is  the  index  of  refraction  structure  constant 
lo  =  inner  scale  of  turbulence 
Lo  =  outer  scale  of  turbulence. 


experimental  verification"that  by  semlemP1ri'ca1  theory  and  direct 


(4)  cn(h)  -  Cn2(H0)  H(n)"4^3 


where  H(n)  is 
The  outer  seal 


path  altitude  divided  by  the  initial  ait-if..w  u 
e  varies  linearly  with  altitude?  aU'tude  Ho- 


(5) 


Lo(n)  =  Lo(Ho>  H(e) 


calculations  allow  their'Sse^or'riarae'Sariet  Tf6  ™ndimenS'onalized 

SS!~  5.SK52  tsF3?—-:' “ 

experimental  quantities.  The  n^  vaJiabljrare6"*  Comb1natio'>s 


l6>  u  ~  *W 


v  =  n/L 

HL  -  Ho 

g  =  -  .  0 


v .  2^0(h0):2 

- c - 

vi  .  2kCL0(HL)J2 


H(n)  =  (1  +  gv) 


simplifie^the^’ntegrals  a^d^phas^zeslh6^  °f  the  path*  This 
outer  scale  dependence  of  theTbulence  statistics"'  ''"1tia1 


(7) 


The  angle  of  arrival  correlation  functions  are  easily  derived  from 
the  phase  structure  function  using  Eqs.  (2a),  (2b),  and  (7). 

C.  Computer  Generated  Curves 

These  calculations  have  been  performed  on  an  IBM  model  360  computer 
using  double  precision  arithmetic,  and  96  point  gaussian  quadrature 
numerical  integration. 

Figure  1  shows  the  normalized  spherical  wave  phase  structure  function 
computed  from  Eq.  (7)  for  an  upward  path  for  two  H^/Hq  ratios. 

The  curves  show  very  little  saturation  as  is  the  case  for  infinite 
outer  scale,  because  the  expanding  upward  beam  sees  an  ever  incr°asing 
L0,  so  that  a  characteristic  beam  width  is  generally  less  than  the  local 
outer  scale.  For  inner  scales  less  than  a  few  centimeters  the  ex¬ 
ponential  term  of  Eq.  (7)  is  essentially  unity  and  therefore  inner 
scale  effects  are  not  included  in  the  phase  structure  function  cal¬ 
culations. 

As  an  example  of  the  use  of  this  figure  consider  a  10p  spherical 
wave  source  propagating  along  a  45°  upward  slant  path  starting  at  an 
elevation  of  1  meter  and  ending  at  an  elevation  of  100  meters  (L  £  140m). 
Further  j^ume  that  the  initial  outer  scale  (L0)  is  2  meters  and  is 

*  oose  p  =  1  meter  .  Hence  the 
)L  is  1.76.  From  the  graph  for 
x  10-3  and  calculate  Ds(l  m)  =  2.11. 
1000  m  (L  =  1414)  then  D_(l  m)  = 

106. 


10-|t*  nr^0  at  one  meter  altitude.  Ch 
normalization  factor  c2(H0)L05/3(H0 
Hl/H0  =  100  and  p/L0  =  1/2  we  get  1.2 
If  the  path  is  lengthened  so  that  Hi  = 


A  similar  computer  calculation  for  a  downward  slant  path  is  shown 
in  Fig.  2,  where  the  curve  has  been  renormalized  to  the  receiver  altitude. 

Figure  3  shows  the  differential  path  contribution  for  H0/H^  =  1000 
for  the  case  considered  in  Fig.  2.  There  we  see  a  large  contribution 
to  the  phase  structure  at  the  lower  altitudes  (end  of  downward  path). 

p 

When  the  calculations  of  Figs.  2  and  3  are  extended  to  large  ^  r^-r 
values,  the  normalized  phase  structure  function  curve  saturates  °v  L 
like  the  horizontal  path  curves  saturate,  and  the  differential  path 
contribution  curve  shows  a  more  rapid  increase  near  the  transmitter  end 
than  at  the  receiver  end  just  as  in  the  horizontal  path  case.  The  data 
for  large  p/L0 (Hl )  are  not  included  because  they  are  beyond  the  range 
of  importance  for  practical  applications. 
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Fig.  2.  Normalized  spherical  wave  phase  structure  function  vs 
normalized  separation  for  a  downward  path. 
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bution  to  spherical  wave 
i  for  a  downward  path. 


Figure  4  illustrates  the  small  aperture  azimuth  angle  of  arrival 
correlation  function9  or  large  aperture  arrival  angle  variance  for 
inclined  and  horizontal  paths  and  Fig.  5  shows  the  same  calculation 
for  a  vertically  downward  path.  The  5.005  km  initial  altitude  is 
greater  than  the  limits  set  by  Reference  6,  but  the  equation  is 
correct  in  the  lower  altitudes  where  the  major  contribution  to  the 
result  is  made  and  it  is  a  reasonable  extrapolation  at  the  higher 
altitudes. 


Fig.  4.  Spherical  wave  small  aperture  azimuth  angle  of  arrival 
correlation  vs  separation  and  large  aperture  mean 
square  angle  of  arrival  vs  aperture  radius. 

Figure  6  shows  the  angle  of  arrival  variance  for  a  slightly  inclined 
path.  The  L~V3  dependence  can  be  determined  analytically  by  dropping  all 
but  the  most  significant  terms  in  the  integral  expression  obtained  by 
using  Eq.  (3)  andfljin  Eq.  (2b). 

D.  Summary 

Typical  examples  of  the  slant  path  spherical  wave  phase  structure 
function,  azimuth  angle  of  arrival  correlation  function  and  angle  of 
arrival  variance  have  been  calculated. 
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SMALL  APERTURE  AZIMUTH  CORRELATION 
LARGE  APERTURE  PHASE  VARIANCE  < 


Fig.  5.  Spherical  wave  small  aperture  azimuth  angle  of  arrival 
correlation  vs  separation  and  large  aperture  mean 
square  angle  of  arrival  vs  aperture  radius  for  vertical 
downward  path. 
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Fig.  6.  Small  aperture  angle  of  arrival  variance  vs  range  for 
an  upward  slant  path. 


27 


E.  Conclusions 


tha+  Ih°  uaria2-0^  the  turbulence  statistics  with  elevation  show 

aJoidinnU?h^f  lifreaJe  in  turbulence  effects  can  be  achieved  by 
avoiding  the  lower  atmosphere  especially  for  the  downward  paths.  y 

-Sr  3r“ ■” ‘"c^'a-rsr^s  -v 

a^ng  sto;5anJSSerf?cafhpa?Ss:eCiPr°Clty  °f  SPhCriCal  wave  Pr0Pa9at1°" 
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in.  OUTER  SCALE  EFFECTS  (N  TURBULENCE  DEGRADED  LIGHT  urAM  SPI  CTRA 
A.  Introduction 

outer  scale  effects><teLoo-amfiturt<>S,'eht,'Um  calcu,atlons  t0  include 

spectra  are  consiJ^ed^o^TSSr^^etrspNe^lLl^a^i:"6"" 

the  flTls'thfaTunis  S*eas?er  r-UV°r  distinct  reasons, 
urements  because  data  need  be  takel  It  ^,strument  ^ime  spectral  meas- 
responding  spatial  corrections ,,r  Jy  °ne  p01nt‘  .The  cor- 
ments  at  many  pairs  of  points  Thp  c  13I  sPec^ra  >^equ i re  measure- 

correlations  and  spectra  IJe  i’n  thLll^  rea!°?  is  that  temporal 
with  theoretical  estimates  to  aain  inf!VeS^SefUi  quantities  to  compan 
lence  parameters.  In  this  context  thpS'Tin*10?  a‘«t  attI,osPhe^'c  turbu- 
averaging  effects  as  do  cnsHaim  they  do  not  suffer  from  aperture 
time  lags  are  large  enough  Jo  JeSrese'nT^n^-^^^6  readil*  available 
than  those  conveniently  available  ThJ,S?h^ai  se|?arations  ™ch  larger 
with  which  to  examine  outer  scaleefSts^^  Pr°Vide  3  9°°d  hand1e 

CHffIfd.r^,?f?^:i9?f^’]^nTatarski.s  (TatarsM,  ,967)  and 

already  given  that  final  1  imiting^xpres^inn^1016^^  similar  to  those 
minimum  of  mathematical  development^  can  be  obtai'ned  with  a 

sphericaJ^waveXcaseCwill,,firstCbeCcont'HnS  J111  be  outlined-  The 
amplitude,  phase  and  phase  considered»  expressions  for  log- 

be  followed  by  simila^derivationffor^he^an61"9  derived-  will 

cases  the  starting  equations  wi?l  hi  I-th  PJ  wave  case-  In  both 
Expressions  for  limitinq TseTllJ  tl T?  ^0"9  W1'th  basic  results, 
the  results  are  discussed  and  conclusioif  drawn  S',bSeqU,!nt  sections 

B-  Calculations 

turbulen«thypothesisaloUthat°whenWe  fI”*  assume  T^'ors  frozen 
v  are  collinear  the  temporal  and  sDatial3  reparation  p  and  wind  velocity 
and  B„s(p) 

r.  * 


0a>  Rx.^)  *  [  h  (Bx,s  (^>t 
for  spherical  waves  and  (Tatarski,  1967) 
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(,b)  Rx,s^  -  Bx,s^> 

for  the  plane  wave  case.  The  subscripts  x  and  s  refer  to  log-amplitude 
and  phase  respectively.  K 

Equations  (la)  and  (lb)  show  how  large  separation  (p)  turbulence 
statistics  may  be  obtained  by  using  the  time  lagged  data  from  a 
moderately  long  record  at  one  point.  Further  the  temporal  phase  dif¬ 
ference  correlation  R6s(t)  is  given  by  v 

(2a)  R6s(t)  =  <Cs(r,t)  -  s(r+p,t)][s(r-VT,t)  -  s(r+p-vx,t)> 


For  a  plane  wave  this  reduces  to 


(2b)  R6s(x)  =  2Bs(vx)  -  Bs(v(t-^)  -  Bs(v(t+^)) 

and  for  a  spherical  wave 

' 

(2C)  |2Bs(iik)-Bs(f(,-f)) 

'  Bs(r(T  +  $r)j’d2 . 

j 

Following  previous  work  the  temporal  spectrum  (Tatarski,  1967)  is 
given  by  ’  ' 

(3)  WXiS(f)  -  4  |"  dt  coS(2,fT)Rx>s(T), 

For  the  spherical  wave  case  the  log-amplitude  (upper  sign) 
and  phase  (lower  sign)  correlation  functions  are  given  by 

(4)  Bx.s<»>  -  BA2f0L  d2  j”d*  «♦„(«)  J0  (¥) 

is.the  1  ight  wavelength.  The  index  spectrum  (Tatarski, 
1961)  to  be  used  is 
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(5) 


Hk)  =  0.033  C 


2  U  1.077  V 

"|U  ) 


+  < 


-11/6 


i?  ^.turbulence  outer  scale.  The  value  of  1.077  for  the 
coefficient  of  1/L0  was  chosen  to  give  two  limiting  forms  for  the 
index  structure  function 

f r2  2/3 
I  cnp 


(6) 


d„(p)  - 


c2  l2/3 


p  «  L, 


p  >>  L 


!"  3  159:log  Plot  of.Dn(p)  the  break  point  determined  by  the 
intersection  of  the  asymptotic  forms  occurs  precisely  at  p  =  L  . 

Following  previous  work  (Clifford,  1971)  Eqs.  (5),  (4),  Hal  are 
substituted  into  Eq.  (3)  and  the  i  integration  is  performed.  The 
derivation  is  continued  by  substituting 

(7a)  z/L  =  (u+1 )/2 

and 

(7b)  a'  =  {(kv/2ttP)2  -  l}e~2 
where 

(7c)  B2  *  1  +  (1.077v/2irfL  )2 

o 

and 

(7d)  a  =  2TrfL1/2/k1/2v 

and  performing  the  a'  integration  in  a  fashion  following  that  of 
Reference  1.  The  result  is 

{8)  Wx,s  =  •132Tr5/2k2/3cn1-7/3v"1fi"8/33"8/3(r(4/3)/r(ll/6) 

x  Re  |  du(  1+  expC-1a2(u2-l)/4]{1F1(l/2;-l/3;ie2n2(u2-l)/4) 

+  (r(-4/3)r(ll/6)/TT1/2r(4/3)(i62n2(u2-l)/4)4/3 


1F1 (ll/6;7/3;  i32fi2(u2-l )/4). 
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(5) 


*(k)  =  0.033  C 


+  <2  J. 


■11/6 


where  Lp  is  the  turbulence  outer  scale.  The  value  of  1.077  for  the 
coefficient  of  1/Lo  was  chosen  to  give  two  limiting  forms  for  the 
index  structure  function 


(6) 


Dn(p) 


r2  2/3 
Cnp 


r2  ,  2/3 

Si  Lo 


«  L. 


p  >>  L 


3  129;log  plot  of.Dn(p)  the  break  point  determined  by  the 
intersection  of  the  asymptotic  forms  occurs  precisely  at  p  =  L  . 

Following  previous  work  (Clifford,  1971)  Eqs.  (5),  (4)  (lal  are 
substituted  into  Eq.  (3)  and  the  t  integration  is  performed!  The 
derivation  is  continued  by  substituting 

(7a)  z/L  =  (u+1 )/2 

and 

(7b)  o'  =  {(kv/2tfP)2  -  1}3~2 
where 

(7c)  B2  =  1  +  (1 .077v/2irfL  )2 

....  0 

and 

(7d)  n  =  2TrfL1/2/k1/2v 

and  performing  the  o'  integration  in  a  fashion  followinq  that  of 
Reference  1.  The  result  is 

(8)  Wx,s  =  •132,r5/2k2/3cnL7/3v"1n"8/3e"8/3(r(4/3)/r(ll/6) 
f1  P  p 

x  Re  j  du(  1+  exp[>ijr(u  -1 )/4]{^F^ (l/2;-l/3;i e2n2(u2-l )/4) 

+  (r(-4/3)r(ll/6)/7T1/2r(4/3)(ie2n2(u2-l)/4)4/3 


1F1 (ll/6;7/3;  i$2n2(u2-l  )/4). 
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where 


-  \  <  arg(x)  <  p 


Equation  (8)  is  given  because  it  is  easily  manipulated  to  give 
asymptotic  forms.  Note  that  the  new  information  is  in  the  factor  a  which 
contains  the  outer  scale.  Indeed  Eq.  (8)  reduces  to  the  corresponding 
equation  in  the  appendixof  (Clifford,  1971)  if  the  outer  scale  becomes 
infinite. 

Using  Eqs.  (9)  to  expand  the  functions  in  Eq.  (8)  in  a  power 
series  in  (u2-l)  and  using  the  stand  J  series  expansion  for  ex,  we 
integrate  with  respect  to  u  to  obtain  the  limiting  forms  for  large 
and  small  n  shown  in  Table  I.  The  expansions  have  been  carried  out 
to  a  number  of  terms  sufficient  to  indicate  the  predominant  behavior. 
Again  if  6  =  •  the  expressions  of  (Clifford,  1971)  result. 

The  sphtrical  phase  difference  spectrum  is  obtained  by  following 
a  similar  procedure.  Inserting  Eqs.  (5),  (4)  and  (2c)  into  Eq.  (3) 
and  performing  the  t  integration  gives 


where  Eqs.  (7)  are  used  to  convert  to  the  new  variables.  With  the 
exception  of  the  6  factor,  this  term  is  identical  with  Eq.  (20)  of 
(Clifford,  1971).  Since  3  ^  1  the  arguments  given  in  that  Reference  for 
dropping  the  second  integral  are  still  valid  with  the  result  that 
W6s(f)  has  the  form  given  in  Table  I  for  large  and  small  f. 


\ 
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TABLE  I 


The  Limiting  Case  Results  for  Plane  and  Spherical 
Wave  Time  Spectra  with  Outer  Scale  as  a  Parameter 


P 


PLANE  WAVE 


log  amplitude 


ns«l  Wx(f) 
f?B»l  W  (f) 

A 

fiB«l  W  (f) 
H3»l  W  (f) 


.850  Cnk2/3L7/3v"1(l-D.45B4/3-1.29B"2/3-.430B‘8/3]n4/3) 
2.19  cjjk^L^VfflB)-873 

phase 

4.38  Cnk2/3L7/3v_1  (OB)-8/3(l-.194(aB)f8/3) 

2.19  c2k2/3L7/V1(ne)_8/3 


phase  difference 


p«  (>L)1/2 


o»  (XL) 


1/2 


w6s(f)  - 


.033  L 
.066  ' 


C2k2Lv5/3 


cos 


^  2irpf 


a 


e-8/3f-8/3 


SPHERICAL  WAVE 


flB«l  Wx(f) 
«3»1  WY(f) 

A 


log  amplitude 

=  .191  C^k2/3L7/3v"1(l-C.646e4/3-.573g"2/3-.19l6“8/3]fi4/3) 
«  2.19  CJkW^3 


phase 

!J3«1  W.(f)  -  4.38  cSkW'V’W-V3  Cl 
ns»l  Ws(f>  -  2.19  cl  k2/3L7/3v'1  (ae)'8/3 


,0436(!!6)8/3: 


p«  (XL)1/2 
P»  (XL)1/2 


W4s'f) 


phase  difference 


r 


/.033l 

[.066J 


C2k2Lv5/3  1 


0-8/3  f~e/3 
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For  the  plane  wave  case,  the  log-amplitude  (upper  sign)  and  phase 
(lower  sign)  correlation  functions  are  given  by 

(ID  Bx>s(p)  =  2A2L  | 


Substituting  Eqs.  (11),  (5)  and  (lb)  into  (3),  performing  the  t  in¬ 
tegration  and  substituting  Eqs.  (7c)  and  (7d)  and 

(12)  z  =  fi2((KV/2irf)2-l) 

gives 


d<  <  JQ(<p ] 


)V  * 


k 

7T 


si 


t) 


Vk) 


w 

(13)  Wx  s(f)  =  .033n2k2/3L7/3C2v"1  |  dz 


.-1/2 


1  +  sin(z+n  ) 

(z+fl2)  J 


(z+A2)' 


A  convenient  closed  form  for  the  z  integration  has  not  been  found, 
+  it  can  be  put  into  a  form  useful  for  limiting  cases  by  employing 


(4) 


sin(x+a)  _ 
(x+a) 


1 


db  cos(b  (x+a)) 


in  Eq.  (13)  and  performing  the  z  integration  to  yield 


W„  ff)  =  .033n5/2k2/3L7/3C^r(4/3)r'1(n/6)!l‘8/3B"8/3. 
x  j  s  n 


(15) 


m2  f  P  / 1.  l.  -l.  2  2^ 

|1F1  ^7*  "  3’  "lb3  n  J 
r(-4/3)r(ll/6)r"1(4/3)1r"1/2(-ibn2e2)4/31F1(ll/6;7/3;-ibeV)j 


1  +  Re  [ 

1  2  ( 
db  eib0  . 

J 

0 

Employing  the  series  forms  for  the  hypergeometric  functions  given  in 
Eqs.  (9)  for  large  and  small  si  then  gives  the  expressions  quoted  in 
Table  I. 


The  plane  wave  phase  difference  spectrum  is  obtained  quite  simply 
by  inserting  Eq.  (2b)  into  Eq.  (3)  and  using  the  shifting  theorem., 
for  a  real  even  function.  The  result  is 


(16)  Wfis(f)  =  2 


1  -  cos 


Ws(f) 


Equation  (16)  is  also  used  in  finding  entries  for  Table  I. 


11/6 
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C.  Discussion 


A  comparison  of  two  limiting  values  for  each  temporal  spectrum  in 
Table  I  reveals  that  only  the  log-amplitude  has  a  significant  change 
in  form  in  going  from  one  limiting  case  to  the  other.  For  the  time 
lagged  phase  and  power  phase  difference  spectra,  the  change  from  one 
limiting  expression  to  the  other  would  probably  be  undetectable  ex¬ 
perimentally. 

The  main  point  of  interest  in  the  inclusion  of  outer  scale  effects 
into  the  temporal  spectrum  calculations  is  the  effect  on  low  frequency 
behavior.  To  this  end  the  low  frequency  limits  associated  with  the 
expressions  given  in  Table  I  have  been  extracted,  and  are  presented 
in  Table  II.  They  apply  for  the  case  where  f  «  1 . 077v/2ttLo,  i.e., 
the  period  of  the  frequency  being  measured  is  greater  than  the  time 
required  for  the  turbulence  field  to  traverse  a  distance  of  one  outer 
scale.  (Extremely  long  averaging  times  or  sophisticated  data  processing 
schemes  may  be  required  in  some  cases  to  obtain  such  data.)  In  any  case, 
an  interesting  point  is  that  the  log-amplitude  and  phase  spectra 
approach  a  constant  limit. 

The  spherical  wave  phase  difference  spectrum  is  of  special  interest 
because  it  is  being  measured  at  RADC.  It  may  be  put  into  a  convenient 
nondimensional  form  by  using  dimensionless  variable  combinations.  The 
form  is 


Figure  1  is  a  log-log  plot  of  Eq.  (17)  when  p  <<  Lo.  The  initial 
slope  is  +2  until  the  1st  break  point  of  =  li^ZZp.,  where  the  slope 
is  -2/3.  The  slope  then  becomes  -8/3  after  the  2nd  Break  point  of 
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41  «  L 
AT  >>  l 

AE  «  L( 
AT  »  L( 

p  «  Al 
p  »  AT 

4[  «  L0 
/AL  >>  L 

o 

4l  «  L0 

4l»  lo 

p  «  at  h 

p  »  AT 


PLANF  WAVE 


TABLE  II 

Low  Frequency  Asymptotic  Results  f  <<  >I‘P77v 

2nLo 


Log  Amplitude 


"x<f>  ■  2-'<A  »■'  (vmf/3 


Phase 


3  Ws(f)  »  4.38C 


tfl 

ws(f)  ■  2.l9cJk2L  ) 


Phase  difference 


«s 


(f }  =  '  8.87}  Cn^Lv'1  ( T°B77)  [l  -  cos(^)] 


SPHERICAL  WAVE 


Log  Amplitude 


/  2irfL  V 

*191 


✓  l  \°' J  r  0 

(f)  .  4.38  CJA,-'  {^m)  ,  . 

o  «  L  l\  ol 

V 17577  ) 


Hs(f)  ■  2.19  CjjlAv"11 


Phase 


Phase  difference 

2ffpf ' 


m 


]») 
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W8s  (f )(»/>) 


D.  Conclusions 


The  previously  published  temporal  spectra  for  infinite  outer  scale 
have  been  directly  extended  to  finite  outer  scales  thereby  verifying 
the  conclusion  that  outer  scale  effects  only  become  important  when  the 
period  of  the  temporal  frequency  under  consideration  is  comparable  with 
or  longer  than  the  time  it  takes  the  wind  to  move  the  turbulence  one 
outer  scale  length.  For  low  temporal  frequencies  an  additional  outer 
scale  dependent  break  point  is  added  to  the  frequency  spectra,  and  hence 
low  frequency  temporal  spectra  can  be  used  to  obtain  the  outer  scale  as 
a  check  of  meteorological  measurements. 
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IV.  SUMMARY  AND  CONCLUSIONS 


To  review,  in  the  final  technical  report  three  areas  were  considered, 
and  one  appendix  was  presented.  The  areas  were  averaging  times  for 
propagation  data,  the  calculation  of  predicted  results  for  various 
slant  path  measurements  and  the  calculation  of  expressions  for  temporal 
spectra  of  several  propagation  functions.  The  appendix  presented  a 
suggested  redefinition  for  the  well  known  coherence  length. 

In  the  area  of  averaging  times  an  examination  of  some  typical 
RADC  propagation  data  showed  that  the  averaging  times  were  sufficiently 
long  to  assure,  that  for  the  mean  square  value,  there  would  be  less 
than  ten  percent  variation  between  ensemble  and  time  averages  and  less 
than  ten  percent  variation  between  any  given  average  and  the  stationary 
value.  The  design  of  a  real-time  averaging-time  meter  to  give  required 
averaging  time  while  data  is  being  recorded  was  also  discussed.  Finally 
the  limitations  to  the  averaging  time  of  temporal  spectra  were  con¬ 
sidered.  It  was  indicated  that  the  scatter  of  spectral  values  is  a 
lesser  restriction  on  averaging  time  than  is  the  frequency  resolution 
desired. 

In  the  area  of  slant  path  calculations,  computer  evaluated  ex¬ 
pressions  for  the  phase  structure  function  and  angle  of  arrival  cor¬ 
relation  functions  for  spherical  waves  were  presented.  These  plots 
show  the  nonreciprocal  effect  of  the  atmosphere  on  beams  propagating 
upward  and  downward  in  the  atmosphere.  Other  curves  of  the  differential 
path  contribution  show  the  exceptionally  deleterious  effect  of  the  lower 
atmospheric  layer  on  downward  traveling  spherical  waves,  indicating  the 
desirability  of  avoiding  when  possible  the  lower  atmospheric  layer. 

In  the  study  of  temporal  spectra,  expressions  given  in  the 
literature  for  spectra  of  log -amplitude,  phase  and  phase  difference  of 
atmospherically  degraded  beams  were  extended  to  include  outer  scale 
effects.  The  results  for  phase  difference  spectra  indicate  for 
example,  a  quadratic  behavior  at  low  frequencies  not  heretofore  known. 

The  appendix  contains  a  recommendation  that  the  coherence  length,  r0, 
be  redefined  in  some  cases  depending  on  the  outer  scale.  Specifically, 
when  the  coherence  length  is  greater  than  twice  the  outer  scale,  it  has 
negligible  effect  on  optical  parameters  and  should  not  be  considered. 

For  values  of  the  coherence  length  less  than  twice  the  outer  scale  it 
should  be  defined  as  before. 
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APPENDIX: 

COHERENCE  LENGTH  REDEFINITION 
A.  Introduction 


coberer?co  length  r0  defined  originally  by  Fried  (Priori  iq^i 

horizontal  oath  of  length  L, Wavelength  7:  2°/l  and  fVT6  °Ver  3 
parameter  c!L  rn  is  qiven  bv  tho  won  ^  and  turbulence  structure 

(6.88/2.91  C^LC§)3/5^  y  the  we  known  exPression  r0  - 

review  concepts  involvinq  r„  and  tn  Tt  !  th?  object  here  to 

definition  to'  sityS  S?J5rTtt. “eflSS  of".  the  , 

might  be  noted.  c“  m  a  t1nite  outer  scale 

length  dell JIl tlon^  1 T’ bl’ ext?"si0"  of  the  coherence 
coherence  length  is  compare  with  or  grlaEef'thaS'Jh^oSt^'scalt 
B.  Review 


defined  and  applied  ^n  thisWWctionWnWW^W  coherence  len9th  was 
assumed.  P  thlS  seLtlon  an  1nfinite  outer  scale  will  be 

Signal-to-Noise  Ratio 

conjunction^with'the'signaUto-noisETratio'of'an  ‘W  l%7)  in 

receiver  detecting  an  atmospherically  degraded ligh^Im 

^ local °oscil''atordbeamatbothh1nci dent  d^ra  ®d  deam  collinear'with 
quantum  efficiency  „  and  ofdiameter  D  ^The^t^  Pho^odetector  of 
for  the  ensemble  average  signal -to-noi^e  ratio  ^/^waVXwnWo^e 
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(1) 


> 


T  =  (e)“j?x  |  E(rl)drl  I 

where 

e  =  electron  charge, 

E(r,)  =  electric  field  of  atmospherically  degraded  beam  at  point  F.  on 
photodetector  surface,  and 

angular  brackets  denote  ensemble  average. 

The  expression  given  in  Eq.  (1)  is  simplified  in  Appendix  A  with 
the  final  result 

(2a)  >  x  J  x  [8  j0  r  dr  T(r)  M(r)j 

where 

r  =  magnitude  of  separation  of  two  points  on  detector  surface 
(2b)  t(r)  =  transfer  function  of  circular  aperture 


(2c)  M(r)  *  atmospheric  transfer  function 
a  exp  (  -  Dw(r) ) 

(2d)  Dw(r)  =  wave  structure  function 
-  Ar5'3. 

The  parameter  A  is  defined  by  Eq.  (2d)  and  will  be  evaluated  subsequently. 

In  the  absence  of  turbulent  fluctuations  or  in  the  limit  of  very 
small  aperture  diameters  M(r)  £  1  and  the  signal -to-noise  ratio  is, 

(see  Appendix  A) 
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A^e„^Slra'-t0-n0iSe  ratf°  value^f  (see  °  =  ' 


1 


(4) 


<S>/N 


i  *  (& 


.  -.88  y' 5 

{ pr) 


s i gna  1  -to-noi strati  o  e^essio^Ej^ta)6  Ixcep^thafth116  und*graded 
diameter,  D  has  been  replaced  by’a  Lco^sLnt^  g?veS  by^"  """ 


(5) 


rQ  =  (6.88/A)3/5, 


inus  r0  is  defined  to  preserve  thp  fS*™  LL£) 

signal -to-noise  ratio,  other  geon,etrieS  are“onsidere°d  InlppST 

Sr^^gHi^at^,5W«B?,sss“ 

we  obtain,  (for  an  infffite  Sier  fcafe)  P  9  ^  ^  a"d  (3) 


(6a)  D 


eff 


r  rD 

8  j  rdr  x(r)  M(r)  |1/2 
r  r  J 


(6b) 


8  [  rdr  - 
)q  " 


cos 


1  fr  \5/3 

7  x  *-**\jr) 


1/2 


terms^of  D^hln^clet  """  -tio,  expressed  in 


(7) 


so  hat  average  signal-to-noise 
aperture  area  *D|ff/4. 


ratio 


is  proportional  to  the  effective 
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scaleDtherrfsreiS?ttl2U^eKdl're-t,i:  ?be  ‘oncept  that  for  Infinite  outer 

S  ^peSe  a,Hr  "  S’9na,-t0-'MiSe  rat1°  ^  -king 

quantities^which'^depend^n'^the^oherence  length. ^  eXpreSS1'ng  other 

Imaging  Resolution 

1magiSe(F?i!lieni966r9trinr?Hj1??  S"fs  1nt0  tbe  theory  of  atmospheric 
throuah  l  i  Co"slder  light  from  a  point  source  propagating 

tnrough  a  turbulent  atmosjihere,  imaged  by  a  lens  of  diameter  n  tn  a  y 

expressed^°^?lP^°^PrI^^S°°r^n°^m^^  ^"^e^ve^e  ol  lW2)  ^ 

he  SDatial  spectrum  is  then  defined  (Fried,  1966)  as  the  resolutioJ,  R. 


(8a)  R  =  jjr(f)  M(f)  df. 


The  resolution  is  more  commonly  written  by  exDressinn  it  in  fDmP 

“  1*?™’  ri  ‘be  input  °r  aperture  pile  ?the  Ze  oJ  tSe™en  ) 

In  Appendix  C  it  is  shown  that  the  spatial  frequency,  7  is  related  to 
the  input  aperture  coordinate  difference,  r,  by7  =  ^-7  and  that  * 
Eq.  (8a)  can  be  expressed  as  *  *  y  a'T*  and  that 


(8b)  R  =  /^~r\D  x(r)  M(r)  rdr 
(Ad.)^  J 0 


1  Tr[)eff 
IT*  5 


(Adi) 

lSeartI°re™syzheCge;;ng?veinrediate,y  ,n  te™ 

(8c)  R  = 

(Ad.) 

Thus  the  resolution  is  also  a  function  of  the  effective  aperture  area. 
Strehl  Imaging  Intensity  Ratio 

at?10SPhen'c  imaging  we  also  consider  the  Strehl  intensity 
atio,  D,  also  called  the  Strehl  definition  (O'Neill,  1963) 

(9a)  D  =  <I(0)>/I  (0), 

centpr’nf ritl0  °f  ,the. ensemble  average  intensity  <I(0)>  in  the 

center  of  til  1? !  6  P°1-t  ™age  to  l°{0)*  the  Intensity  in  the 
center  of  the  point  image  in  the  absence  of  atmospheric  degradation. 
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Equation  (9a)  is  expressed  in  terms 
terms  of  aperture  plane  coordinates 
familiar  form 


of  image  plane  coordinates, 
(see  Appendix  C)  it  takes  on 


In 

a 


(9b)  D  = 


in  terms  of  the  effective  aperture  diameter  this  becomes 
/  D  \2 

(9c)  0  =  (-SJL.)  . 


fora  JfetherrelolI!ti<!n’tR  lIldYl50  I-™5  °Ut  t0  be  a  "°™alited 
diameter.  t10"’  R-  a"d  a  function  of  the  effective  aperture 

Imaged  Spot  Size 

imaging,  the  imagedbspJtas?zeer(Bracke^y  1™'?"'’“?^°"  "1th  atmsplieric 
terms  of  the  ensemble average’, • ntensit^’  Th!t.,1s  defi'  ed  in 

•P?»4,r«,.zrit£  dtes?1 A  M  rr- %  "spot 

inVF?g9 /f?  ^hus^eTve5  ^  e"Semble  intensi^  IslhL 

J* 

(10a) 


^o 

~r 


* 

<I(r2)>dr, 

=  <T((T)~ 


Appendix*?™  °f  the  aPGrtUre  Plane  Variables  this  becomes  (see 


16 

Tx 


(Ad.)‘ 


dOb) 

77  8jM(r)t(r)r  dr 
and  finally  in  terms  of  the  effective  diameter, 

Ooc)  d  >  fidL 
0  . 
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Focused  Spot  Size 

.  .  iFin?11?  °neuCan  consider  the  size  of  a  spot  focussed  through  the 

turbulent  atmosphere  onto  a  receiving  plane.  The  focused  spot  size 

&w!!ibe  dfinJdT^ctly  as  the  imaged  spot  size  in  terms  of  an 
integral  scale  of  I(p2),  the  intensity  in  the  receiving  plane. 

o  r  r 


(11a) 


Trdj:  ||l(p2)dF, 

T"=  ~TToT 


It  is  shown  in  Appendix  D  that  the  expression  for  I(p0)  is  formally 

Dlane1diJ^nth  f°r  fhe  tmage  plane  intensity  whe6  the  image 
plane  distance  d-j  is  replaced  by  the  range,  L.  The  resultinq  ex¬ 
pression  for  the  focussed  spot  size  is  then  9  ex 


(lib) 


iLhqr1i'ffcat?LfmC-SHd/Ht0  im??ad  Sp0t  sizes  df/do  <s  then  the  geometrical 
magmrication  m  -  df/d0  -  L/d^  as  one  would  expect. 


in  terms  of  the^ffective  anerture'dia^t™"51?^611  Can  be  exPressed 

s  arssaSig 11 -ssio" 

»■*«  -  »ter 

C.  Examination  of  D  « 

- _ cff 

iMidfSrtBSMSriSTS^S  tL'cKJg  &£;  T  *5 

t'zszr&i  s^s^a-s:  ™r  ™V 

of  JWf  FS  slmiffSS! I,!J!5  1,E?  cote™  length,  the  behavior 
*p.3lrf»  “  'id,  JiJrsi:  «K°;sh  ;9  th4t  of  th,:  s™n 

Thf 

be  thIhsameeflsni?n°Eqsf  SCale  effects  1,111  *>™.lly 

extended  to  include  outer  scaJe'effects.  WaV<!  StrUCture  fu"ct™ 


02)  0iff  *  8 


rdr  — 


•  ?  Dw<r> 


spectrum  whlcMncfude^  a  turbulence  spatial 


(13)  *(K)  = 


0.033  C‘ 


(oeJ  •  ■’) 


Mn)d=fi^^3,WandCfo0re?  *>  L^DJr)™  V*  LQ. 

antes  because  of  the  second  cr?ter?in!  as  3i5cussed1n  Append?!  ^;°77 
waves^re^ive^b^Eqs^tlAa^and6, 14b) ^respectively6  ^  spher1cal 
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(14a)  D  ( r )  =  871^  L  [1  -  J_(k  r)]$  (k)k  d< 
wp  j  q  U  II 

L  2  °° 

(14b)  Dws(r)  -  8,2k2|o  d,  g  }o  Cl-J0(Kr)>(^)  K  dK 


These  two  functions  are  plotted  in  Fig.  2.  They  have  the  expected 
behavior.  For  small  values  of  r 

040  Dws(r)  =  |  Dwp(r)  =  |  2.41  k2L  C2  r5/3 

Both  of  these  will  be  represented  by  the  infinite  outer  scale  ex¬ 
pression 

<14d)  Dwp,s  ■  6.88(r/r0)5/3  -  6.88(r/rop>5)5/3 

where  r0p»  and  ros  are  given  in  Eqs . (B7 )  and  (BIO).  For  large  values 
of  r  they  both  saturate  at  the  same  value  which  is  shown  in  Appendix 
E  to  be 

04e)  Dws>pH  -  1.381  k2LC2  if3 

The  fact  that  the  wave  variance,  Bw(0)  =  ^  Dw(°°),  was  the  same  f°r 

both  plane  and  spherical  documents  the  fact  that  the  mean  square  log 
amplitude  plus  phase  fluctuation  generated  along  a  single  line,  con¬ 
necting  transmitter  and  receiver  depends  only  upon  the  path  and  not 
upon  the  type  of  wave  either  side  of  the  path. 

We  can  now  also  write  the  infinite  separation  value  Dw(°°)  in 
terms  of  the  coherence  length  for  plane  and  spherical  waves,  re¬ 
spectively. 

/L  \5/3 

(15a)  DwpH  -  3.265^J 

/L  \5/3 

05b)  DWSH  =8.71  f^J 
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The  limiting  expressions  for  Dw  will  now  be  used  in  the  calculation 
of  asymptotic  forms  for  the  effective  diameter. 


To  provide  a  feeling  for  how  the  effective  diameter  should  behave 
tor  various  ranges  of  the  actual  aperture  diameter  we  qualitatively 
examine  the  integral  expression,  Eq.  (11),  for  the  various  cases  of 
interest.  The  procedure  will  be  to  sketch  curves  for  the  various 
actors  in  the  integral  definition  to  demonstrate  the  limiting  cases 
The  two^factors  of  interest  are  the  optical  transfer  function, 

T(r)  =  “  cos  (fi\-  -  (r/D)2  which  goes  from  1  at  r  =  0  to  0  at 

•  v  •  ■  --------  /  1 

exr>  j  -  j 

ffor  r  <<  L  and  assumes  the 


r  =  D  and  the  atmospheric  transfer  function  M(r) 

(  /  \  5/3) 

which  drops  off  as  exp  -  /r  '  1 

limiting  value  of  exp 


h  m 


Vr>) 


{ 


3. 


for  r  >>  L 


The  value  of  the 


,  §  < 

constant  C  depends  on  the  type  of  wave  as  reflected  in  Eqs.  (15). 
Sketches  of  -r(r)  and  M(r)  are  shown  in  Fig.  3a.  From  Eqs.  (12), 
the  effective  diameter  is  then  the  square  root  of  area  under  the  curved 
formed  from  r  times  the  product  of  the  two  curves. 


We  now  sketch  the  two  curves  for  various  ranges  of  D. 

Case  I.  D  «  rQ  <<  LQ  (see  Fig.  3b). 

This  is  the  infinite  outer  scale  small  aperture  case  considered  in  the 
previous  section.  M(r)  is  essentially  constant  for  0  <  r  <  D  so  that 


(16a) 


fD 

'0 


dr  r  T(r)J^ 


=  D. 


Case  II.  rQ  <<  D  <<  Lq  (see  Fig.  3c). 

This  is  the  large  aperture  case  considered  earlier.  T(r)  is  essentially 
constant  over  the  range  where  M(r)  is  significant,  so  that 


(16b) 


rdr  exp  {-3.44(r/rQ)5//3 


r 


o 


Case  III.  D  >>  rQ  >>  L0  (see  Fig.  3d). 


This  case  is  new.  In  this  case  M(r)  drops  to 
before  r(r)  goes  to  zero  so  that 


an  asymptotic  value 
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(17a)  M(r)  %  exp  {-  7  C(LQ/ro) 5/"3  } 
over  almost  all  of  the  range  0  <_  r  <_  D  and 

(17b)  Deff  %  [exp{-  ^  C(L0/rQ)5/3}  x  8  |°  drr  i(r)]1/2 

=  D  exp{-  }  C(L0/r0)5/3}  =  D  exp{-  }  Dw  (-)}  . 

Thus  Deff  is  again  proportional  to  D,  this  time  the  constant  of  pro¬ 
portionality  being  exp(-  1/4  Dw  (“)).  These  three  regions  are  shown 
in  the  leg-log  sketch  of  D  ff  vs  D  in  Fig.  4. 


The  significant  point  here  is  in  the  behavior  in  region  III 
For  small  values  of  the  ratio  U/r0  the  sidestep  as  determined  from 
the  factor  exp  {-  C/4  (LQ/r0)573}  ;s  negligible  and  the  line  in 
region  III  is  a  continuation  of  the  line  in  region  I.  For  that  case 
the  behavior  is  no  different  from  the  infinite  outer  scale  case 
discussed  in  the  first  section.  However,  for  larger  values  of  Ln/rn, 
the  line  in  region  III  is  sidestepped  significantly  and  the  behavior 
in  the  two  regions  is  different. 


.  .  P16  d,|^re|r|ce  in  the  two  regions  is  indicated  by  the  amount  of 
sidestep  obtained  from  the  expression  exp  {-  C/4  (Ln/rn)5/3}  This 
quantity  is  plotted  as  a  function  of  L0/r0  in  Fig.  5  for  both  plane  and 
spherical  waves.  There  we  see  that  for  L0  equal  to  r0  the  sidestep  is 
1/0.4  =  2.5  for  plane  waves.  (Actually  the  factor  exp  {-  C/4  (L  /r  )5/3} 
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gives  a  change  in  the  vertical  position  which  is  then  translated  into  a 
sidestep  of  1/exp  {-  C/4  (L0/r0)5/3}.  The  corresponding  value  for 
spherical  waves  is  1/0.10  for  spherical  waves.  For  the  outer  scale  one 
half  the  coherence  length,  the  straight  line  portions  in  regions  I  and 
III  of  the  Deff  vs  D  graph  should  be  within  a  factor  of  2.5  of  each  other 
for  both  plane  and  spherical  waves.  This  will  be  used  in  the  next  section 
where  a  suggested  extension  of  the  definition  of  the  outer  scale  will  be 
presented. 


D.  Extension  of  the  Definition  of  r 
_  o 


In  the  previous  section  we  examined  the  limiting  cases  of  the 

behavior  of  the  effective  diameter  as  a  function  of  outer  scale  and 

coherence  length,  and  noted  that  with  finite  but  large  outer  scale, 

Deff  saturates  at  the  value  r0  as  it  did  with  infinite  outer  scale. 

Indeed  with  outer  scale  less  than  half  the  coherence  length,  the  curve 

of  Deff  versus  D  would  show  less  than  a  factor  of  two  variation  in  going 

from  D  less  than  Lq  to  D  greater  than  Lq.  This  behavior  is  indicated 

more  clearly  in  Fig.  6  where  we  see  computed  plots  of  Deff  vs  D 

covering  the  range  of  D  much  less  than  L  to  much  greater  than  L 

o  o 
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EFPECTIVE  DIAMETER  O^tt 


APERTURE  DIAMETER 


Fig.  6.  Effective  aperture  diameter  Deff  vs  actual  diameter 
for  various  values  of  ratio  outer  scale  to 
coherence  length  (L  /r  )*  for  a  spherical  wave. 
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for  several  ratios  of  L0/rp.  We  see  that  for  L0/r0  <  1  the  transition 
from  one  straight  line  portion  of  the  Deff  versus  D  curve  to  the  other 
is  nearly  undetectable.  Therefore  a  representative  generalization 
when  L0/r0  £  1  is  that  the  Deff  versus  D  curve  is  essentially  the  no 
turbulence  case,  where  Deff  =  D. 

Thus  it  is  recommended  that  the  coherence  length  not  be  defined 
for  cases  of  outer  scale  less  than  one  half  the  coherence  length,  and 
that  for  cases  with  outer  scale  greater  than  one  half  the  coherence"" 
length  the  definition  of  coherence  length  for  infinite  outer  scalp  hp 
retained. This  extension  is  the  formal  embodiment  of  the  fact  that  the 
coherence  length  is  a  useful  quantity  for  r0  <  L0  but  that  for  r0>  Ln 
the  behavior  is  essentially  that  of  the  no-turbulence  case. 

E.  Discussion 


There  are  a  few  loose  ends  yet  to  pick  up.  The  first  is 


-  -  -  -  «««  up.  me  1 1  rb l  i s  a  dis¬ 

cussion  of  the  physical  behavior  for  the  r0  >  L0  case.  Basically  the 
meaning  is  contained  in  the  wave  variance  Bw(0)  =  .69  LC^  i_5/3  which 
is  the  maximum  value  of  Dw(r)/2.  In  terms  of  Bw(0)  the  larges?  ratio 

°en^rD  1/2  Bw(0)).  When  the  maximum  phase  variance  is  less 

than  0.5  rad  ,  i.e.,  0.707  rad  RMS  deviation  or  less  than  a  quarter  wave, 
then  the  optical  effects  of  the  atmosphere  are  essentially  negligible 
In  such  a  case  the  resolution  will  be  unimparied.  For  such  a  case  a 
quantity  which  measures  the  atmospheric  resolution  degradation  will 
cease  to  have  itJ  full  meaning  and  can  be  safely  neglected. 

F.  Summary  and  Conclusions 


In  this  appendix  we  have  considered  the  effects  of  atmospheric 
turbulence  on  optical  heterodyne  detection  and  on  optical  resolution 
and  spot  size.  Specifically  we  considered  the  case  when  the  outer 
scale  of  turbulence  is  comparable  with  the  coherence  length,  r0,  and 
suggested  an  extension  of  the  definition  of  the  coherence  length. 

First  the  heterodyne  signal-to-noise  ratio  and  the  imaging 
quantities  were  shown  to  all  be  functions  of  the  effective  aperture 
diameter,  a  quantity  which,  for  infinite  outer  scale,  is  equal  to 
the  aperture  diameter  for  very  small  apertures,  and  is  equal  to  the 
coherence  length  for  very  large  diameters. 


The  effective  aperture  diameter  was  then  investigated  for  finite 
outer  scales.  It  was  found  that  when  the  outer  scale  was  small  enough 
to  be  less  than  one  half  the  coherence  length  that  the  effective  aper- 
turG  diameter  varied  with  actual  diameter  essentially  as  if  there  were 
no  turbulence.  Thus  it  is  recommended  that  when  the  outer  scale  of 
turbulence  was  less  than  half  the  coherence  length  that  the  coherence 
length  be  undefined  and  that  when  the  outer  scale  is  larger  than  the 
coherence  length  that  the  usual  definition  of  coherence  length  applyl 


56 


REFERENCES 

Brackey,  T.A.,  (1968),  "Preliminary  Report  on  the  Use  of  Spatial 

lD0VflJrifyi968°PTheaOh?^rfCa1V-nS  ^st0"s'”  Report  2156-9, 
Department  of ’tlortH^iSJat-  Un”;ers’tJ'  ElectroScience  Laboratory, 

BP  N  NO. "l&SST  2d 

Fo^Mel  0MoF°rCe  SySt™S  Comnand'  ^Hht-Patterson  Air 

,rt“r  fciTwSF £s  «“•“"*  &— 

°,Nc1ikie;Lp;b{!S  nptirs- Add1s°"- 


57 


APPENDIX  A 


Eq  m  tI1FnaPmdl’X^arC  pres?nted  the  steps  necessary  to  reduce 
tq.  (I)  to  Eq.  (2)  and  to  examine  certain  limitina  cesec 
the  ensemble  average  square  of  the  intearal  in  f<-9  n f!;  W,r  1  9 
integral  and  exchanging  integrations  and  summation  gives'  *  ^  6 

(Al)  <||E(r1)d?1  |2>  =  <  jj  E(rj)  E*(Fy)  > 

aperture 


IS 


<E(rj)E*(r1")>drjdr^'  . 
aperture  * 

^tch  sum  and  difference  coordinates  r  =  FT  r»  and  ^  = 

Iri  +  ri)/2-  Usi"9  the  def1n1t’°"  of  the  atmospheric  mutual  coherence 
function,  (Fried,  1967), 

(A2)  <E(r-])E*(r-j')>  -<|E|2>  M(|Fj  .  p.|)  =  M(FK 

and  using  the  aperture  function  f (Fj ) ,  given  by 
(A3)  f(r,)4,:  |F,'-D 

°i  !r] I  >  D  , 

we  have 


m  <l/E<’-i)d'-1|2>=|//|f(Fj)f(P1')<E(Fj)E*(P)')>  dF.’dP; 

-.00 

00 

jj  drh(r)  jj  dR  f(1T  +  £  )  f(  IT  -  £)1. 

-°°  L  -oo  J 


(A5) 


<  I E  1 2> 


T(r)  byG  1nner  int69ra1  15  related  t0  the  aPerture  transfer  functi 


(A6) 


on 


oo 

II  dir  f(K  +  £  )  f  (  IT  -  r  )  =  T(r)  fj  f(f)  =  nD^  T(r) 

-00  j  j  A 
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For  a  circular  aperture 
(A7)  <(r)=f[cos-1$-S 

Inserting  Eq.  (A6)  into  Eq.  (A5)  gives 

2°° 

(A8)  <|  E(r, )dr, |2>  =  <|E|2>  ^2-  dF  M(r)T(r) 

1  ^ 

mOO 

Since  M(r)  and  t(r)  are  functions  of  r,  magnitude  of  F,  we  switch  to 
polar  coordinates  and  perform  the  angular  integration,  giving 

(A9)  <|  E(r-|)dr.|  |2>  =  <|E|2>  r  dr  M(r)  x(r). 

Equation  (A9)  is  the  desired  result.  Inserting  Eq.  (A9)  into 
Eq.  (1 )  gives  Eq.  (lb)! 

Equation  (A9)  can  be  examined  for  two  limiting  cases.  For 
aperture  diameters  small  compared  with  the  coherence  length  M(r)  is 
constant  over  the  region  of  interest  and 


This  is  used  in  the  derivation  of  Eq.  (2). 

For  aperture  diameter,  D,  much  larger  than  the  coherence  length, 
i(r)  =  1  over  the  region  of  interest,  and  Eq.  (A9)  becomes 


(A12) 


Ef^  )dr-|  |2>  = 


( A 1 3 ) 


3  r (6/5) 
5 
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Equation  (A14)  is  used  in  the  derivation  of  Eq.  (4). 
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APPENDIX  B 

In  this  appendix  we  present  general  expressions  for  coherence 
length  for  plane  and  spherical  waves  for  cases  where  is  a  function 
of  range.  First  consider  the  plane  wave  case 

r\  r\  f  L  f°° 

(Bl)  D  (r)  =  87TZkZ  dz  d<  <*(<)  [1  -  Jn  (<p)] 
w  j0  JQ  n  o 

Taking 

(B2)  $(k)  =  0.33CZ(2)k_1  1/3 

and  making  the  substitutions 
(B3)  u  =  kp 

and 

(B4)  v  =  z/L 

gives 

1  oo 

(B5)  D  (r)  =  .033  x  8A2Lr5/3  dv  C2(vL)  du  u"8/5  [1  -  J  (u)] 
w  jQ  n  jQ  o 

Further,  from  Eqs.  (2d)  and  (5)  of  the  main  text. 

(B6)  Dw(r)  =  6.88  ^  3 

Using  known  values  for  the  u  integral,  (Tatarski,  1961,  p.  269)  and 
solving  for  rQ  gives  the  desired  results 
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For  a  spherical  wave 


(B8) 


Dw(p)  =  8 


A2I>(^!>  -  j0<«)] 


substituting  for  4>(k)  from  Eq.  (B2),  for  <  and  z  from  (B3)  and  (B4) 
respectively,  and  simplifying  gives 


(B9) 


]  °° 

Dw(p)  =  0.033  x  8A2Lp5/3  dv  v5/3  C2(Lv)  du[l-J  (u):u"8/3 

Jq  n  Jn  0 


Evaluating  the  u  integral  and  solving  for  rQ  using  Eq.  (B6)  gives 

l3/5 

(BlOa) 


ro 


.2  . 


6.88 

2.91  k2L 

Tdu  ub/3  C2(Lu) 

and  when  Cn  is  a  constant 


(BlOb) 


ro  = 


6.88 


13/5 


(2.91 )*(3/8)  k2L  C2  , 
n  J 


Equations  (B7)  and  (BIO)  are  the  last  results  of  this  appendix,  general 

formulas  for  the  coherence  length, (r0) when  is  a  function  of  longitudinal 
position.  /m 
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APPENDIX  C 


In  this  appendix  we  present  derivations  of  several  expressions  for 
imaging  situations.  The  first  is  the  derivation  relating  the  normalized 
Fourier  transformation  of  the  image  plane  intensity  I(r2;  to  the 
spectrum  i (f)M(f)  and  also  to  the  integral  over  the  input  aperture  in 
Eg.  (6b).  The  spectrum  of  the  image  plane  intensity  is  given  by 


(Cl) 


P(f)  = 


ii27rT-?2  <I(7?)>  df2  . 


In  Eq.  (Cl)  7  is  spatial  frequency  in  cycles  per  meter. 

The  image  plane  intensity  is  related  to  the  input  plane  fields  by  the 
well  known  relationship 

•  k  -  - 


(C2) 


(C3) 


<I(r2)>  =(^)  <  IE(rl}  6 


HT  1  2 


dr 


1 


‘(A} 


dr. 


1  5-  (r,'-r.") 


I.l,l  1 

dr1,,<E(r1,)E+(r1,,)>  e  1 


Following  steps  identical  to  those  in  Appendix  A  we  switch  to  the  same 
sum  and  difference  coordinates,  express  <E(r-| 1  )E(r-|  ")>  in  terms  of  the 
atmospheric  transfer  function,  and  perform  the  integration  over  the 
sum  coordinate,  ]f.  The  result  is  comparable  to  Eq.  (A8)  and  reduces 
to  Eq.  (A8)  for  r9  =  0.  The  result  is 


(C4) 


<I(r2)>  = 


=te) 


<|E| 


ttD 

~T 


dr  T(r)  M(r)  e 


k 

d. 


inserting  this  expression  for  I(r2)  into  Eq.  (Cl)  and  performing  the 
r2  integration  gives 


(C5) 


dr  T(r)M(r)6 


We  note  that  P(0)  can  be  easily  evaluated  since  M(0)  =  t(0)  =  1.  The 
result  is 


(C6) 


P(0)  =  < | E | 2> 


T“ 
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»'"Lthe  r  i"te9ra‘fon  <"  Eq.  <C5)  and  normalizing  the  result 
to  Waives  one  desired  result,  the  form  of  the  normalized  triform 


£{5}  ■  r(>d.f)  Hlxd^f). 


The  integral  expression  for  the  resolution 
by  integrating  Eq.  (C5)  with  respect  to  f 
in  the  delta  function.  The  result  for 


can  be  quite  eas_ily  obtained 
and  noting  that  f  occurs  only 


» 


PT^f 


is 


(C7) 

(C8) 


R=f^dT=(:~)  |  dr  T(r)M(r) 

2tt  fO 

=  7—2  r  dr  M(r) 

(Ad.)4"  Jq 


Eq.  (C8)  is  another  desired  result. 


To  obtain  the  aperture  plane  expr 
ratio,  we  use  Eq.  (C4)  and  simplify  by 
tegration.  The  result  is 


ssion  for  the  Strehl 
performing  the  angul 


intensity 
ar  r  in- 


(C9) 

It 

O 

t— 1 

II 

O 

2tt 

\  rdr  x(r)M(r) 

yoi 

2tt 

rdr  T(r) 

To  obtain  the  imaged  spot  diameter  as 
use  Eq.  (C4).  Thus  since  t(0)  =  M{0)  =  1 


defined  in  Eq.  (10a)  we 
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(Cl  2) 


fej 


2  2  d2 

<IEI  >  V"  x  2n  I  rdr  T(r)M(r) 


Eq!eO0b5  the'dis'?lda?esin!)  1nt°  Eq-  (,°a)  and  reamm91 


ng  gives 
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APPENDIX  D 


In  this  appendix  we  consider  a  light  beam  that  is  focused  through 
a  turbulent  atmosphere  to  a  beam  waist.  An  expression  is  derived  for 
the  focused  spot  size  in  that  situation.  The  development  parallels 
that  of  deWolf  (deWolf,  1971)  but  the  definition  is  one  which  reduces 
to  proper  values  in  the  case  of  vanishing  turbulence  and  finite  aperture. 


Let  the  transmitting  aperture  and  observation  planes  be  designated 
respectively  by  transverse  coordinates  p-j  and  p2  and  longitudinal  values 
z 1  and  zo.  Then  assuming  the  limitations  of  Fresnel  diffraction,  Eo(p9) 
the  field  at  the  observation  plane  is  given  by  i 


(D-l) 


E2(p2) 


ik  e 


2 
P? 

-ik(L  +  ^  ) 


2ttL 


JJ 


dpj  E(p.j)  exD 


where 


exp 


(=* 


(p^-2p.j  *p2 


B( 


Pi  >p J 


(D-lb) 


E(p-|  )  exp(-ikp^/2L) 


is  the  transmitting  aperture  plane  field 


(D-lc)  L  =  z2  -  z1 


(D-lc)  B ( p-j  ,p2)  =  modification  to  Greens  function  due  to  turbulent 

atmosphere. 

The  ensemble  average  irradiance  pattern,  I(p2)  in  the  focal  plane,  R=L 


(D-2)  I(p2)  =  <|E2(F2)|2>  - 

\2 


2  *j 

(?x)  SSII  dpldpl  E^1  )E*(P]'  >  eL  p2  *  (pTpP  <B(p1,^2)B*(pj,p2)> 


Assuming  an  isotropic  reciprocal  atmosphere  (DeWolf,  1971), 

<B(pi,p2)  B*(pj,p2)>  is  the  spherical  wave  atmospheric  transfer 
function 

(D-3)  <B(prp2)  B*(F1',p2)>  =  M(|7j-  P2|)  =  exp(-^Dws(  Ip^-p’-j  |) 
where  Dws(  |p-]-pj  I)  is  the  spherical  wave  structure  function. 
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Writing  Eq.  (D3)  in  terms  of  sum  and  difference  coordinates 
(D-4a)  r  =  ^ 


(D-4b) 


R  = 


+  i' 


gives 


(DS)  I(p2) 


(^l)  I 


ik 


dr  e 


-^Dws<r>+T 


dR 


E1(^+  £  )  E^(R-  £  ) 


The  inner  integral  is  related  to  the  aperture  transfer  function, 
t(F).  For  the  usual  case  of  a  roundoaperture  of  diameter  D  illuminated 


by  light  of  constant  intensity, 
Eq.  (A-5)  by 


E-|  |2,  the  inner  integral  is  given  in 


( D— 6 ) 


where  -r(r)  is  the  aperture  transfer  function 


(D-7) 


and 
(D— 8) 


I(p2) 


dr  M(r)  x(r)  e 


i  k  p 
r  p2 


r 


The  focused  spot  size  is  defined  by  the  normalized  integral  scale 
as  was  the  imaged  focused  spot  size  (see  Appendix  A) 


(D-9) 


~r 


1  ( p  2 )  d  p  2 

"TO) 


This  definition  has  the  advantage  that  it  is  applicable  to  gaussian 
beams  as  well  as  evenly  illuminated  apertures,  and  reduces  to  the 
standard  finite  well  definied  value  for  an  evenly  illuminated  aperture 
in  the  absence  of  turbulence. 
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fn^_,f?ua^on.  !D"?)  is  simply  evaluated  by  noting  that 

d  andY  Eq‘  (C”4)  if  we  interchange  ~2 

i  nJ  ih?4  the  ™me™tor  and  denominator  of  Eq7 
been  evaluated  in  Eqs.  (C-10)  and  (C-12)  respectively 
expression  for  Eq.  (D-9  is  peLiiveiy. 


Eq.  jD-8)  is 
and  r2  and 
(D-9)  have 
The  resulting 


(D-10) 


4aL 

ttTT 


eff 


the  focused  spot  size  is  m,  where  p  1 


m  = 


Vdi 

37-r- 


# 


m  is  merely  the  geometrical  magnification  as  one  might  expect. 
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APPENDIX  E 


In  this  appendix  we  present  the  integrals  which  determine  the 
numerical  constants  in  the  Von  Karmann  index  spectrum.  Let  us  assume 
that  the  spectrum  has  the  form 


(E-l)  =  C2n  A/((B/L0)2  VV1/6 


k  is  the  magnitude  of  the  three-dimensional  spatial  frequency. 
We  determine  the  constants  A  and  B  by  the  requirements  placed  on  the 
associated  structure  function,  Dn(r) 

(E-2)  Dn(r)  =  C*  r2/3  r  «  Lq 

(E-3)  Dn(r)  -  C l  l*f3  r  »  LQ 


This  choice  is  made  so  that  on  a  log-log  plot  of  Dn(r)  vs  r,then 
extensions  of  these  asymptotic  lines  intersect  at  r=L0  giving  the 
effective  breakpoint  at  that  value. 

For  an  isotropic  spectrum,  the  spectrum  and  structure  function  are 
related  by 

oo 

(E-4)  Dn(r)  =  8„  ^  A  - 

Inserting  Eq.  (E-l)  into  Eq.  (E-4)  and  making  approximation  (E-2) 
gives 


(E-5)  cnV'3  =  8,  c3  A  jo  (i  -  iis-srj  k-5/3  dK. 

Using  (Tatarski,  1961,  p.  270),  we  have 

(E— 6)  C2  r2/3  =  8tt  C2  Ar2/3  x  7r/2r  (8/3)  sin  (tt/3) 

where  r(x)  is  the  gamma  function.  Thus  we  have 
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(E-3)  giv2stin9  Eq'  (E_1)  int°  Eq*  (E"4)  and  makin9  “PP^wliMtlons 


K2dK 


(E_8)  Cn  Lo/3  =  cn  8"  a  f  - - ... 

'0  «B/L*)V)11/6  * 

Using  (Gradshteyn,  1965,  Eq.  3.251-2) 

(E-9)  =  C„2  8,A  IT2'3  L2/3  \  B  (f  ,  ’  ).  C2  8,A  B‘2/3L2/3  1  ill 

V  '  0  r(ll/6) 

where  B(x,y)  is  the  Beta  function.  Using  the  identity  r(z)r(l-z)  = 
tt  esc  irz ,  we  find  J  '  '  v  ' 


(E-10)  B 


.{5  r(i)]3 

J 


=  1.077. 


The  final  result  is 


(E-ll)  $n(K)  =  0-033  Cn 


iwl 


+  K 2  J 


\n/6 


Now  consider  the  plane  and  spherical  wave  structure  functions  for 
finite  outer  scale  and  large  separation.  From  Eq.  (B-l)  we  have  for 
a  plane  wave  '  r 

00 

(E"12)  Dw,P£(p)  =  j  d*  I  d<  K$n(K)  [1  -  Jo(Kp)] 


J0  '0 

and  from  Eq.  (B-8)  for  a  spherical  wave 


(E_13)  Dw,sp(p)  =  87l2|<2 


lofT) 


1  -  Jq(kp) 


Substituting  <'  =  K  in  Eq. 
letting  p  be  large  so  that 
both  equations,  it  is 


(E-12)  and  <*  =  kL/z  in  Eq.  (E-13)  and 
Uq (kp )  <<  1  gives  the  same  expression  for 
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(E-14)  Dw(«)  =  87i2k2 


dz 


d<‘  k'  (k1) 


Substituting  for  <tn(ic)  from  Eq.  (E-ll)  and  performing  the  integrals 
gives  for  both  plane  and  spherical  waves 


(E-15) 


DW<")  ■ 


x  0.033  x  ( VS)k2 


(1.077)5/3 


L  C *  L5/3 
n  o 


=  1.381. 


This  value  is  quoted  in  Eq.  (14e)  of  the  text. 
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